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SUMMARY 


r 

For  a  nurr.ber  of  turbulent  flow  problems,  a  first-order  closure 
or  eddy  transport  approach  to  the  modeling  of  the  turbulence  is 
not  sufficiently  accurate.  It  is  necessary  that  a  more  powerful 
approach  be  used  to  study  the  dynamics  of  the  turbulence.  Flows 
involving  chemical  reactions,  such  as  those  in  a  chemical  laser, 
are  good  examples  of  flows  requiring  a  more  powerful  approach.  The 
mixing  in  HF  and  DF  chemical  lasers  takes  place  under  conditions  of 
very  large  heat  release.  The  coupling  between  the  heat  release  and 
the  turbulence  is  an  important  feature  of  the  flow.  A.R.A.P.  has 
been  funded) AFO-SR  Contracts  F4 4620-73-U-0027  and  P  4 1*6-2 0-75-C- 
QQ2f*to  employ  this  basic  method  to  investigate  problems  connected 
with  chemical  lasers  of  particular  interest  to  the  Air  Force  Weapons 
Laboratory.  These  studies  have  resulted  in  considerable  insight 
into  the  complex  processes  occurring  in  laser  cavities. \  These 
studies  have  been  documented  in  References  1,  2,  and  3- 

”  'This  s  report  describes  recent  results  of  studies  using 

the  reacting  shear  layer  (RSL)  computer  program  and  discusses  the 
proposed  procedure  for  handling  multi-step  chemical  reactions  in  a 
second-order  closure  computation  of  turbulent  reacting  flows. ^ 
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I .  INTRODUCTION 


A  second-order  closure  program  for  the  investigation  of  tur¬ 
bulent  ,  compressible,  chemically  reacting  flowfields  Is  under 
development  at  A.R.A.P.  Second-order  closure  methods  have  provided 
an  understanding  of  the  physics  of  turbulent  flows  which  could  not 
be  provided  by  classical  first-order  closure  techniques  such  as 
the  eddy  viscosity  method.  We  have  successfully  applied  these 
techniques  to  a  variety  of  turbulent  flow  problems  including  incom- 
pressible  and  compressible  shear  layers/  the  behavior  of  the 

ll 

planetary  boundary  layers,  and  the  dispersal  of  chemically  reacting 

g 

pollutants  in  a  turbulent  atmosphere.' 

A.R.A.P.  is  funded  under  AFOSR  Contract  No.  FM620-75-C-0026 
to  develop  an  invariant  second-order  closure  model  capable  of  investi¬ 
gating  the  mixing  and  reactions  in  HF  and  DF  chemical  lasers.  The 
mixing  in  these  devices  takes  place  under  conditions  of  large  local 
heat  release  due  to  the  highly  exothermic  chemical  reaction.  The 
interaction  between  the  heat  release  and  the  turbulence  is  an  impor¬ 
tant  feature  of  the  flow. 

The  set  of  second-order  closure  equations  governing  the  mixing 
and  reaction  of  compressible ,  turbulent,  multi-component  reacting 
shear  layers  has  been  modeled  and  programmed  for  solution  in  our 
RSL  program.  This  interim  scientific  report  describes  some  recent 
results  of  studies  using  this  program  and  discusses  the  proposed 
procedure  for  extending  the  program  for  studying  multi-step  chemical 
reactions  in  a  second-order  closure  computation. 


II.  REACTING  SHEAR  LAYER  (RSI.)  PROGRAM 


The  reacting  shear  layer  program  has  been  completely  checked 

out  over  the  past  six  months.  Program  predictions  for  nonreacting 
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flows  5  (two-dimensional  wake,  heated  wake,  heated  planar  jet) 
are  in  good  agreement  with  experimental  measurements.  Comparison 
of  program  predictions  with  experimental  data  for  reacting  flows 
is  now  underway.  The  RSL  program  has  been  used  for  computations 
of  the  HP  chemical  laser  flowfield,  hydrogen-air  premixed  combustion 
and  diffusion  flames  and  propane-air  premixed  combustion  and  dif¬ 
fusion  flames.  The  calculations  have  been  performed  using  two  dif¬ 
ferent  models  for  the  scalar  correlations.  One  model  (used  during 
the  model  development  stage)  was  to  set  all  third-order  (and  higher- 
order)  scalar  correlations  to  zero.  The  other  model  is  the  A.R.A.P. 

’’typical  eddy”  box  model.  Preliminary  results  for  exothermic  react- 
7 

ing  flows  suggest  that  there  Is  only  a  small  difference  between  the 
calculations  using  these  two  models.  More  extensive  tests  are 
necessary  to  check  the  validity  of  this  result  for  different  flow 
geometries  and  in  different  flowfield  regions.  It  is,  of  course, 
obvious  that  if  the  model  involving  setting  the  higher-order 
correlations  to  zero  can  be  used  for  most  of  the  calculations,  it 
v/ill  lead  to  considerable  program  simplification  and  increase  in 
speed . 

The  use  of  the  ’’typical  eddy”  box  model  in  these  calculations 
has  verified  the  basic  concept  of  modeling  the  scalar  correlations 
using  a  probability  density  function  composed  of  a  set  of  Dirac 
functions  of  variable  strengths  and  position  in  the  scalar  phase 
space.  The  computations  have  also  allowed  us  to  determine  some 
shortcomings  of  the  present  model.  A  number  of  simple  modifications 
have  already  been  made  to  the  model  to  extend  its  applicability  to 
more  severe  environments ,  and  some  other  changes  are  currently  under 
study . 


The  current  version  of  the  program  is  restricted  to  one-step 
forward  reactions.  In  the  past  month,  we  have  developed  a  procedure 
for  handling  multi-step  reaction  systems  including  the  complete 
turbulence-chemistry  interaction.  The  procedure  is  briefly  described 
in  the  next  section. 


III.  SECOND-ORDER  CLOSURE  FOR  MULTI-STEP 
REACTION  SYSTEMS 

The  "typical  eddy"  box  model  has  been  most  extensively  used 
up  to  now  with  three  species  (a,  B  and  y)  flowfields.  In  principle, 
the  model  is  quite  general  and  can  be  extended  to  a  large  number 
of  species.  However,  the  direct  addition  of  each  new  species  in 
che  second-order  closure  system  will  add  ten  equations  for  the 
means  and  second-order  correlations  and  does  not  appear  to  be  a 
reasonable  procedure  for  complex  reaction  systems.  Actually,  such 
a  procedure  is  not  necessary. 

Consider  the  mixing  and  reaction  of  two  initially  separate 
streams  o  and  B  .  These  streams  can  be  known  mixtures  of  a 
large  number  of  chemical  species.  The  streams  mix  and  react 
following  a  large  number  of  elementary  reactions  and  form  a  mixture 
of  product  species  designated  as  y  .  We  assume  that: 

1.  The  composition  of  mixtures  a  and  6  is  fixed  for 
the  entire  flowfield.  There  are  no  internal 
reactions  among  the  species  in  these  two  reactant 
mixtures . 

2.  All  the  product  species  that  compose  y  are 
molecularly  mixed.  There  are  internal  reactions 
within  the  y  species  mixture.  Due  to  the  dif¬ 
ferent  rates  of  the  various  elementary  reactions 
(between  species  present  in  a  and  3  mixtures 
and  the  species  present  in  y),the  composition  of 
the  mixture  y  varies  at  different  points  across 
the  flowfield. 

Consider  the  following  two-reaction  system  as  an  illustration 
of  the  procedure.  The  complete  set  of  reactions  that  we  will  be 
using  in  our  study  of  the  DP  chemical  laser  is  presented  in  Table  I. 
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P0  +  F  -i-  !)!•’  +  P 

i\  +  n  pf  +  v 

a.  P_,  known 

B  He  +  F  +  F.,  known  proportion:' 

il 

Y  He  +  P.,  +  F,  +  P  +  F  +  PF 


i 


3  (He  +  F  +  F„) 


y(He  +  P-  +  F_  +  P  +  F  +  DF) 


cx  ( D  )  — 


The  interaction  between  turbulence  and  chemistry  only  has 
to  be  taken  into  account  For  reactions  between  a  and  3  , 
a  and  y  ,  and  3  and  y  .  Species  within  y  are  assumed  to 
be  molecularly  mixed  and,  therefore,  species  correlations  do  not 
have  to  be  considered  for  the  internal  reactions.  For  conven¬ 
ience,  we  rewrite  the  reaction  system  as  riven  below,  tarninr 
each  chemical  species  with  the  mixture  that  it  is  a  part  of: 

D.,(a)  +  F  ( 3 )  PF(y )  +  D(y) 

el 

D,(a)  +  F(y)  +  PF(y)  +  P(y) 

c 

D2(y)  +  F(  3 )  -  PF(y )  +  D(y) 


(y )  +  F ( y ) 


PF(y)  +  P(y)  internal  reaction 


F;,(fl)  +  P(y)  -  V'F(y)  +  F(\) 

F,(y)  +  IH  Y )  -*•  PF(y)  +  F(y)  internal  reaction 


Wo  can  now  write  the  chen’.i.'ul  source  tenaa  for  the  component 
species.  For  example, 


dn.,(rt ) 

dt 


P . ,  ( a )  F  (  R ) 
M«)F(y) 


+  P;(a)F’(P) 


+  n;(rt)F’(Y) 


dP0(y) 

dt 


IMy)P(S)  +  ld(Y)F'(a) 

ITT  (  y  )  f  (  y  ) 


d  C  . 


Two  other  pieces  of  information  are  now  inquired  to  complete 
this  formulation.  First,  we  will  he  solving  equal  lens  for  the 
second-order  eorrela  t;  1  ons  ct 1 3 f  ,  otryr  ,  and  F’y1  .  Now, 

D,’ (a)P’  (B)  ,  Pi  ( a ) F r ( y ) >  etc.,  can  be  simply  calculated  from 

f-  L. 

these  correlations  when  one  knows  the  composition  of  the  mixtures, 
and  so  the  source  terms  are  completely  determined.  Second,  wo  have 
to  find  the  total  loss  of  5  and  P  .  In  the  above  example,  f 
is  being  lost  due  to  the  consumption  of  both  F(6)  and 
We  determine  the  term  that  causes  the  greater  loss  of  P  and  use 
the  figure.  If  at  some  point  the  consumption  of  F(P)  leads  to 
the  larger  loss  of  P  ,  then  one  has  to  keep  track  of  the  corres¬ 
ponding  components  F0(P)  and  He( P)  that  now  become  part  of  y  . 
The  composition  of  Y  has  to  be  calculated  at  each  point  across 
the  flow field . 

This  procedure  operates  within  the  framework  of  three  over¬ 
all  species  mixtures  a  ,  P  ,  and  Y  and  the  models  and  equations 
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